Bosonization in higher dimensions via noncommutative field theory 
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We propose the bosonization of a many-body fermion theory in D spatial dimensions through a 
noncommutative field theory on a (2D — l)-dimensional space. This theory leads to a chiral current 
algebra over the noncommutative space and reproduces the correct perturbative Hilbert space and 
excitation energies for the fermions. The validity of the method is demonstrated by bosonizing a 
two-dimensional gas of fermions in a harmonic trap. 

PACS-98: 04.70.Dy, 11.15.-q 



Bosonization of one-dimensional fermion systems is a 
well-established technique in condensed matter and field 
theory [1,2]. Its usefulness lies in the fact that the low- 
energy dynamics of the equivalent bosonic theory encodes 
collective excitations of the fermion system and thus pro- 
vides a handle to analyzing strongly correlated fermions. 

The success of the method relies on special proper- 
ties of one-dimensional spaces. Intuitively, the fermion 
dynamics around the Fermi (or Dirac) sea become tan- 
tamount to one-dimensional wave propagation, the cor- 
responding phonon states encoding quantum excitations 
[3] . This leads to a mapping of states between the two 
systems both at the many-body and the field theory (op- 
erator) level. 

There are several approaches to extending these tech- 
niques to higher dimensions [4,5]. The most obvious way 
to proceed is by generalizing the Fermi sea idea 'point- 
wise' in the extra dimensions, giving rise to the so-called 
radial bosonization. This approach, however, has some 
obvious drawbacks that are at the heart of the difficul- 
ties with higher dimensions. Specifically, it produces one 
tower of states for each point in the extra dimensions, 
thus overcounting the degrees of freedom by not taking 
into account the 'finite cell' structure of the sea. Improve- 
ments to the method have been proposed, but bosoniza- 
tion of higher dimensional fermion systems remains to a 
large extent an open issue. 

The purpose of this Letter is to propose a theory that 
addresses the above problem and provides a potentially 
exact bosonization scheme for higher dimensional theo- 
ries. The approach is through a phase space description 
of the many-body fermion problem. For D spatial dimen- 
sions this leads to a 2D — 1-dimensional nonlocal bosonic 
theory. (In this respect, our approach is closest to the 
one in [5].) Only in one dimension the dimensionalities 
agree and the theory becomes local. We accept this as 
the price of the (nontrivial) payout of higher-dimensional 



bosonization. 

We shall present here an outline of the new method, 
motivated by the semiclassical phase space droplet pic- 
ture, and demonstrate it in a simple, solvable situation. 
This would make its physical meaning and content as ac- 
cessible as possible. A full derivation will be given in a 
subsequent publication. 

We shall consider a collection of spinless, noninter act- 
ing fermions in D = d + 1 spatial dimensions. (Spin 
and interactions can be added later without conceptual 
problems.) The starting point will be the semiclassical 
phase space droplet description of the system [6] , in the 
hamiltonian formulation derived and fully exposed in [7] . 
We review here the basic notions and quote the relevant 
results. 

Particles move on a phase space of dimension 2d + 2 
with coordinates <^> M and a single-particle Poisson struc- 
ture 



{^,r}s P = e fW , n,v = l,...2d + 2 



(1) 



For simplicity we shall assume that we have chosed Dar- 
boux coordinates and momenta 4>^ — {x l ,p 1 } so that B^ v 
is in the canonical form, with determinant det d^ v = 1. 
Particle dynamics are fixed by a single-particle hamilto- 
nian H sp ((f)^). 

The system can be described semiclassically by a state 
in which each cell of volume (2irTi) d+1 in phase space 
is filled by a single fermion. This leads to a phase space 
'droplet' picture of constant density p a = l/(27r?i) d+1 , de- 
termined by its 2d + 1-dimensional boundary. We choose 
to parametrize the boundary in terms of one of the phase 
space coordinates, say p 1 , which becomes a function of 
the remaining variables R(<t>^) (0 M ^ p 1 )- It is conve- 
nient to separate in notation the variable conjugate to 
p , x 1 = a ({p , x } S p — — 1); and use middle Greek in- 
dices for the remaining 2c? phase space variables 4> a and 
corresponding Poisson structure af3 . cf> will collectively 
denote the coordinates d> a . 



1 



The dynamics of the system are given in terms of the 
canonical Poisson brackets of the field R(u, 4>) (not to be 
confused with the single-particle brackets {., .} sp ). Using 
the shorthand Ri = R(<Ji, <pi), these are derived as [7]: 

{R 1 ,R 2 } = — [-5'(a-)6(4>-) -Slp-^daR+dfiStt-)] 

Po 

(2) 

with (cr_,0_) and (a + ,(f> + ) relative and mid-point coor- 
dinates respectively. The hamiltonian for the field R is 
the sum of the single-particle hamiltonians over the bulk 
of the droplet 

H = Po J dp l dad 2d cf)H sp {<j, <P)§{R - p 1 ) (3) 

where i?(a;) = |[1 + sgn(x)] is the step function. (2) and 
(3) define a bosonic field theory (in a hamiltonian setting) 
that describes the system semiclassically. 

The above theory arlready represents an improvement 
over previous methods, which essentially amounted to 
keeping only the first (constant) term in the Poisson 
brackets. Althought at this stage it still suffers from 
shortcomings similar to the ones of previous Fermi sea 
approaches, it will serve as the starting point for an exact 
bosonization. Below we shall simply motivate the form 
of the exact theory and demonstrate its merits, leaving a 
full derivation and details for a future publication. 

We start by noticing that R+ in the right hand side 
of the above Poisson brackets can also be put to Ri (or 
R2), since the antisymmetry of 9 al3 makes the difference 
irrelevant. In view of this, the right hand side can be 
recast using the single-particle Poisson brackets for the 
variables <fif as 

^daR+dpStt-) - {R{a u <Pi), <5(0i - 4> 2 )} sp i (4) 

We now make the central claim of this paper: the 
exact theory will be obtained from the above expres- 
sion by 'quantizing' (at the classical level) the remain- 
ing phase space coordinates; that is, by promoting 
the single-particle Poisson brackets to noncommutative 
Moyal brackets {.,.}* on the 2d-dimcnsional phase space 
manifold </>": 

{R U R 2 } = — [-<5'(<r_)5(0_)-<y(<r_){i2i,<y(0_)}.i] 

Po 

(5) 

The Moyal brackets between two functions of (f> are 
expressed in terms of the noncommutative Groenewald- 
Moyal star-product on the phase space 4> [8] : 



Gm* = ^ \F{<t>) * G(4>) G{<j>) * F{4>)] (6) 

with h itself being the noncommutativity parameter. 
Correspondingly, the hamiltonian H is given by expres- 
sion (3) but with star-products replacing the usual prod- 
ucts between its terms. 

Inherent in the definition of the star product is a one- 
to-one mapping between classical functions F(<p) and op- 
erators Fona first-quantized phase space. Specifically, 
define the set of quantum operators <f> corresponding to 
the residual classical phase space, satisfying 

[$ a ,^] =ih6 afi , a, (3=1,..., 2d (7) 

and choose a basis \a) for the (single-particle) quantum 
mechanical Hilbert space on which the above operators 
act. This could be, e.g., the Fock space of d harmonic 
oscillators, or any other basis. Classical commutative 
functions map to operators by adopting a specific opera- 
tor ordering : F(4>) : in the expression for F(<j>), and star 
product maps to operator product: 

F(<f>) <-» F = : F{4>) : , (F * G)(<f>) <-»: FG : (8) 

It is customary to choose the fully symmetric Weyl or- 
dering 

. e i,"*" := e UfV , qeR 2d (Q) 

in which case the corresponding star product is given in 
terms of the Fourier transform of the functions as 

{p * dm = J (0^ (fc " ^^ h ^ qakf> ( 10 ) 

Any other appropriate ordering can be used, leading to 
alternative definitions of the star product. 

In the limit h — > the above becomes the usual (com- 
mutative) convolution integral in Fourier space, while 
the brackets (6) go over to standard Poisson brackets on 
the manifold. For nonzero h the star-product and cor- 
responding Moyal brackets are nonlocal. The obtained 
theory is, therefore, nonlocal in the 2c? dimensional phase 
space while it remains local in the a direction. 

To analyze the content of the above theory, it is bene- 
ficial to perform a change of variable in the field R to a 
matrix representation. We define the mapping between 
functions F(<f>) and matrix elements F ab = (a\F\b): 

F ab = f d 2 UF^)C ab {4>) , F{4>) =J2 Fa " C ba(<P) 

J a,b 

(11) 
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with C a b(4>) an appropriate set of basis functions depend- 
ing on the ordering: 



C ab (<P) = 



d 2d q 



"(a\ :e- iq * : |6) (12) 



In this representation, star products become matrix mul- 
tiplication and phase space integrals become traces: 

(F * G) ab = F ac G cb , j d 2d <t>F(<f>) = (2nh) d F aa (13) 

with repeated indices summed. 

The above mapping can be applied to the classical dy- 
namical field R(a, 4>) mapping it to dynamical matrix el- 
ements R ab (a). To express the full Poisson brackets for 
R in this representation we need the matrix expression 
for the delta function 5(0i — 4*2), with defining property 



d^0iF(^i)5(0! 



F(<h) 



(14) 



Since <5(</>i — ^2) is a function of two variables, its matrix 
transform in each of them will produce a symbol with 
four indices S aibl '' a2b ' 2 . The above defining relation in the 
matrix representation becomes 



(2Trh) d F aibl 5 biai;a2b2 — F a2b2 
which implies 



caibi;a 2 b2 



<ra 1 b 2 ra 2 bi 



(15) 



(16) 



With the above, and using p = l/(27r7i) d+1 , the 
canonical Poisson brackets of the matrix R ab become 

{Rf ,Rf } = - 2Trh5'(a_)6 ad 5 cb 

+ 2m6(a- ) {Rf5 cb - R cb 5 ad ) (17) 

Upon quantization of the theory, the fields R ab {a) be- 
come operators whose quantum commutator is given by 
the above Poisson brackets times ih. Defining, further, 
the Fourier modes 



Rf = J 



da 
2^h 



R ab (<j)e 



-ika 



(18) 



the quantum commutators become 

[Rf, R cd ] = kS(k + k')8 ad 5 cb - Rf +k ,5 cb + Rf +k ,5 ad 

(19) 

The Casimir i?g a = N is the total fermion number. For 
a compact dimension a, normalized to a circle of length 
2ir, the Fourier modes become discrete. 

This is nothing but a chiral current algebra for the 
matrix field R^ on the unitary group of transformations 



of the first-quantized states \a). To make this explicit, 
consider for the moment that the first-quantized Hilbcrt 
space is if -dimensional, that is, a,b,c,d — 1, ... K (this 
would be the case for a compact phase space {<p a }), 
and define the hermitian K x K fundamental genera- 
tors of U(K), T A , A = 0, ...K 2 - 1, normalized as 
tr{T A T B ) = \b AB and defining the U(M) structure con- 
stants [T A , T B ] = if ABC T c with f AB = 0. Using the 
T as a basis we express the quantum commutators (19) 
in terms of the R A = tr(T A R) as 

[Rt R*\ = \m + k')5 AB + l f ABC R c k+k , (20) 

This is the so-called Kac-Moody algebra for the group 
U{K). 

The coefficient k KM of the central extension of the Kac- 
Moody algebra (the first, affine term) must be quantized 
to an integer to have unitary representations. Interest- 
ingly, this coefficient in the above commutators emerges 
quantized to the value k KM = 1. This is crucial for 
bosonization [2]. The k KM — 1 algebra has a unique irre- 
ducible unitary representation over each 'vacuum'; that 
is, over highest weight states annihilated by all R A (k) 
for k > and transforming under a fully antisymmet- 
ric SU(K) representations under T A (0). We argue that 
these Fock-like representations correspond exactly to the 
perturbative Hilbert space of excitations of the many-body 
fermionic system over the full set of possible Fermi sea 
ground states. The U(l) charge which is a Casimir, 
corresponds to the total fermion number; diagonal oper- 
ators R k , for k < and H in the Cartan subgroup of 
U{K) generate 'radial' excitations in the Fermi sea; while 
off-diagonal operators R k , for k < and T off the Car- 
tan subgroup, generate transitions of fermions between 
different points of the Fermi sea. 

Finally, the hamiltonian of the bosonic theory becomes 



H 



^L tr H sp (<j,p\mR-p 1 ) (21) 



where p 1 remains a scalar integration parameter while <j> 
and R arc matrices as before. Clearly there are issues of 
ordering in the above expression, both quantum and ma- 
trix, just as in standard 1 + 1-dimensional bosonization. 

It is clear that for D = 1 our theory reproduces stan- 
dard bosonization. In this case there is no residual phase 
space and we recognize R(cr) as the derivative of the 
bosonized field d a <& representing the chiral density in one 
dimension. (Two independent chiral fields R± , represent- 
ing the two intercepts of the compact droplet boundary 
with the line x=constant, are needed.) The expression 
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of the hamiltonian (21) for H sp = p or H sp = p 2 /2 as- 
sumes the usual relativistic or nonrelativistic bosonized 
form H = (d x <S>) 2 /2 or H = (d x <P) 3 /6, respectively, fully 
recovering the standard results. 

To demonstrate the applicability of this theory we shall 
work out explicitly the simplest nontrivial example of 
higher-dimensional bosonization: a system of N nonin- 
teracting two-dimensional fermions in a harmonic oscil- 
lator potential. The single-particle hamiltonian is 

H sp = ^{p 2 + x 2 +p 2 +x 2 ) (22) 

where for simplicity we chose the oscillator to be isotropic 
and of unit frequency. The single-body spectrum is 
the direct sum of two simple harmonic oscillator spec- 
tra, E mn = h(m + n + 1), m,n = 0,1,.... Calling 
m + n + 1 = K, the energy levels are Ek = hK with 
degeneracy K. 

The TV-body ground state consists of fermions filling 
states Ek up to a Fermi level Ep = UK p. In general, 
this state is degenerate, since the last energy level of de- 
generacy Kp is not fully occupied. Specifically, for a 
number of fermions N satisfying 

N = Kf( - K £ Z 1 ) + M , < M < K F (23) 

the Fermi sea consists of a fully filled bulk (the first term 
above) and M fermions on the Kp-degenerate level at 
the surface. The degeneracy of this many-body state is 

g{K F ,M) = M ^p- M)\ < M > 

representing the ways to distribute the M last fermions 
over K F states, and its energy is 

E(Kp,M) = n K ^ KF -^ 2KF - 1 KnKpM (25) 

6 

Clearly the vacua (Kp, M = K F ) and (Kp + 1, M = 0) 
are identical. Excitations over the Fermi sea come with 
energies in integer multiples of h and degeneracies ac- 
cording to the possible fermion arrangements. 

For the bosonized system we choose polar phase space 
variables, 

1 x- 

hi = -{ Pi + xf) , 0i = arctan — , i = 1, 2 (26) 

2 Pi 

in terms of which the single-particle hamiltonian and 
Poisson structure is 

{9i, hj} sp — Sij , H sp = hi + /12 (27) 

For the droplet description we take hi — R and 9i = a 
which leaves {/i2,#2} ~ {^2,^2} as the residual phase 
space. The bosonic hamiltonian is 



H = j^-^ J dadh 2 d6 2 (^R 2 + h 2 R) (28) 

The ground state is a configuration with R + h 2 = Ep. 
The nonperturbative constraints R > 0, h 2 > mean 
that the range of h 2 is < h 2 < Ep. 

We 'quantize' the single-particle residual phase space 
(h 2l 2 ) by defining oscillator states |a), a = 0,1,2,... 
satisfying h 2 \a) = h(a+^)\a). The nonperturbative con- 
straint for h 2 is implemented by restricting to the Kp- 
dimensional Hilbcrt space spanned by a — 0,1,... Kp 
with Ep = hKp — 1. In the matrix representation 
R ab becomes a U(Kp) current algebra. We also Fourier 
transform in a as in (18) into discrete modes i?° b , n = 
0, ±1, . . . (a has a period 27r). The hamiltonian (28) has 
no matrix ordering ambiguities (being quadratic in R and 
h 2 ) but it needs quantum ordering. Just as in the 1 + 1- 
dimensional case, we normal order by pulling negative 
modes N to the left of positive ones. The result is 

j = Yl R -n R n + \« + (a + \)K a (29) 

n<0 

To analyze the spectrum of (29) we perform the change 
of variables 

Rn = K-a+b + (a-K F + l)S ab S n (30) 

The new fields R satisfy the same Kac-Moody algebra as 
R. The hamiltonian (29) becomes 

TT 1 1 

\ pot iyba 1 1 not fyba , / tv- \ jjaa 

-J-- 2^ H -n H n + 2 H H + ( K F - ^)n 

n<0 

| Kp(K F - l)(2K F - 1) (3i) 

The above is the standard quadratic form in R plus a 
constant and a term proportional to the U(l) charge 
R% a = N - Kp{K F - l)/2. 

The ground state consists of the vacuum multiplct 
\Kp,M), annihilated by all positive modes R n and 
transforming in the M-fold fully antisymmetric irrep of 
SU(Kp) (0 < M < K F - 1), with degeneracy equal to 
the dimension of this representation Kp\/M\{Kp — M)\. 
The U(\) charge of R is given by the number of boxes in 
the Young tableau of the irreps, so it is M. The fermion 
number is, then, N = K F (K F - l)/2 + M. Overall, we 
have a full correspondence with the many-body fermion 
ground states found before; the state M — Kp is ab- 
sent, consistently with the fact that the corresponding 
many-body state is the state M = for a shifted Kp. 
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The energy of the above states consists of a constant 
plus a dynamical contribution from the zero mode Rq. 
The quadratic part contributes \UM, while the linear 
part contributes Ti(Kp — \)M . Overall, the energy is 
TiK F {K F - 1){2K F - l)/6 + KK F M, also in agreement 
with the many-body result. 

Excited states are obtained by acting with creation 
operators R- n on the vacuum, introducing n quanta of 
energy. Due to the presence of zero-norm states, the 
corresponding Fock representation truncates in just the 
right way to reproduce the states of second-quantized 
fermions with an SU(Kf) internal symmetry and fixed 
total fcrmion number. These are, again, into one-to-one 
correspondence with the excitation states of the many- 
body system, built as towers of one-dimensional excited 
Fermi seas over single-particle states -E m>ra , one tower for 
each value of n, with the correct excitation energy. We 
have the nonperturbative constraint < n < Kp, as well 
as constraints related to the non-depletion of the Fermi 
sea for each value of n, just as in the one-dimensional 
case. The off-diagonal operators R!^ create particle tran- 
sitions between towers, with the total particle number 
fixed to N by the value of the Rq. 

We conclude by stating that the proposed bosonized 
theory is defined in a hamiltonian setting. The transition 
to a lagrangian requires 'inverting' the canonical Poisson 
brackets to a two-form u> and corresponding connection 
uj = dA. The resulting action will contain the noncom- 
mutative generalization of the Wess-Zumino term, ex- 



pressed in terms of a star-unitary field U given by the 
star-exponential of a real field U = exp„(i$) [9]. This 
will be presented in a future publication. 
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